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Abstract. The nonlinear interaction, due to quantum electrodynamical (QED) 
effects, between two electromagnetic pulses and a radiation gas is investigated, ft 
is found that the governing equations admit both modulational and filamentational 
instabilities. The instability growth rates are derived, and the results are discussed. 
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Within quantum electrodynamics (QED), there are new and interesting phenomena with 
no classical counterparts. A prominent example is the Casimir effect. Similarly, photon- 
photon scattering (see, e.g., [1-4]) mediated by virtual electron-positron pairs, does no 
occur within classical electrodynamics, in which electromagnetic waves in vacuum are 
indifferent to each other. The collisions of photons with photons have attracted much 
interest over the years, both from an experimental and an astrophysical point of view 
(see [5-20] and references therein). The effect of photon-photon scattering could be of 
fundamental importance in high-intensity laser pulses, in ultra-strong cavity fields, in 
the surrounding of neutron stars and magnetars, and in the early Universe. 

In this paper we shall investigate the interaction between two electromagnetic pulses 
via a background radiation gas. The gas acts as a mediator of the interaction, and the 
pulses exchange their energy through this background. Moreover, in contrast to the 
case with no radiation background, even parallel propagation of the pulses results in 
non-zero coupling. It is furthermore shown that the resulting system of equations lead 
to both modulational and filamentational instabilities. Here we present the growth rates 
for these processes. Our results are finally discussed in the context of astrophysical and 
laboratory applications. 

We express the nonlinear self-interaction of the photons in terms of the Heisenberg- 
Euler Lagrangian as found in Ref. [1] 

L = e F + ml [AF 2 + 7G 2 } , (1) 

where F = (E 2 — c 2 B 2 )/2 and G = cE ■ B. Here re = 2a 2 h 3 /45m 4 e c 5 ss 1.63 x 
10 _30 ms 2 /kg, a is the fine-structure constant, h is the Planck constant divided by 2n, 
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m e is the electron mass, and c is the speed of light in vacuum. The Lagrangian (JTJ) is 
valid when there is no electron-positron pair creation, and when the field strength is 
smaller than the critical field, i.e. 

uj <C m e c 2 /h, \E\ <C E crit = m e c 2 /e\ c (2) 

respectively. Here e is the elementary charge, A c is the Compton wavelength, and 
E crit ^ 10 18 V/m. 

Following Ref. [5], we write the dispersion relation for a low energy photon in a 
background electromagnetic field (E, B) as (see also Ref. [15] and references therein) 

u{k, E, B) = c\k\ (l - |A|Q| 2 ) , (3) 



where 



\Q\ 2 = £o 



E 2 + c 2 B 2 — (k ■ E) 2 — c 2 (k ■ B) 2 — 2ck ■ (E x B) , (4) 



and A = A±, where A + = 14k and A_ = 8k for the two different polarisation states of 
the photon. Furthermore, k = k/k. The approximation A|Q| 2 <C 1 has been used. 

We will below study (a) plane electromagnetic pulses propagating in a background 
consisting of a radiation gas in equilibrium, and (b) a radiation gas affected by an 
electromagnetic (EM) pulse propagating through the gas. 

The electromagnetic field in Eq. (@J) is here a superposition of the fields due to the 
two coherent electromagnetic pulses, with electric fields Ei^ir, t) exp[z(fei j2 -r — o;i j2 t)] + 
c.c., and the incoherent field due to the photon gas. Here c.c. denotes complex conjugate. 
The amplitudes of the pulses are supposed to be slowly varying as compared to their 
rapidly oscillating phases. In the absence of the pulse, the photon gas is assumed to 
be isotropic. The validity of the Lagrangian requires that the electromagnetic fields 
are smaller than the critical field. Thus we can use a perturbative approach where the 
photon gas is only weakly perturbed, i.e., to lowest order in a wave amplitude expansion, 
the photon gas is still isotropic. 

Now we need to establish the connection between the dispersion relation ([3~|) 
(originally derived for static or slowly varying background fields), and the contribution 
from the dynamical radiation fields of our problem. Here it is crucial that the QED 
nonlinearities do not involve the derivatives of the electromagnetic fields. In particular, 
we may write the general wave equation resulting from the Heisenberg-Euler Lagrangian 
as 

(4^- v2 )^ = J2 [kkiE 3 E k E l + g llkl B'B k E l + h ijkl B^E k E l ) .(5) 

The specific forms of the tensors /, g and h can be found by comparison with the 
wave equations in for example Refs. [12, 14]. We note that the nonlinear terms are 
cubic, i.e. the lowest order scattering processes involve four waves. The coefficient 
\Q\ 2 in (@J) results from carrying out all possible summations, where one of the fields 
in the right hand side is taken as the wave field, and the others as the constant 
fields, and finally taking the scalar product with the polarization vector. However, 
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for combinatorial reasons, the coefficient of interaction becomes the same for waves 
influenced by background fields that are dynamical waves. This can be seen by noting 
that the combination surviving the wave mixing is obtained by taking one field as the 
original wave, the other as a wave field at a different frequency and the third as its 
complex conjugate t.. Other combinations surviving the wave mixing are ruled out by 
the conservation relations. 

Next we substitute the particular fields into (j3J) in order to evaluate their effects on 
the dispersion relation (JHJ). Here we note that the rapidly oscillating contributions to \Q\ 2 
should not be included, i.e., terms proportional to E 1 -E 2 exp[i((k 1 + k 2 )-r — (ui+u 2 )t)}, 
E\ exp[2i(fex -r — Uxt)) etc. can be neglected, since in Eq. (JHJ) such terms will be averaged 
out. Thus the only terms kept in the expression for \Q\ 2 are proportional either to the 
photon gas density or to the slowly varying intensity in one of the pulses. We note that 
our calculation of the effect of the pulses on the photon gas requires the intensities to 
vary on a long scale as compared to the wavelengths of the gas photons. 

For case (a), the relations (fcj ■ E) 2 = E 2 /3, (k t ■ B) 2 = B 2 /3, and E x B = hold, 
where fcj, i — 1,2, stands for the pulse wavevectors. Hence, from (J3J) we obtain 

\Q 1 \ 2 =\gg + ax£ 2 , (6) 

|Q 2 | 2 = |4 + «2^i, (7) 

where S g = e (E 2 + c 2 B 2 )/2 and = e \Ei\ 2 is the energy density of the radiation gas 
and pulse i, respectively. Using ej = Ei/\Ei\, we have defined 

a 1<2 = 2 - 2k x ■ k 2 - (ki,2 • e 2 ,i) 2 - [(^1,2 x k 2 ,i) ■ e 2>1 )] 2 . (8) 

In case (b), the directions k g of the photons in the gas are approximately random. 
The last statement holds true as long as the pulse energy densities are small, thus 
inducing only weak deviations from the equilibrium values of the radiation gas quantities. 
The error in this approximation is therefore only of first order in the pulse energy 
densities, giving a second order correction in the final equations. Then (JJJ) yields 

IQ/ = !(<£ + <&+/V^), (9) 

where we have introduced the parameter 

(3 = ei • e 2 + (fei • k 2 )ei ■ e 2 - • e 2 )(k 2 ■ e x ). (10) 

With the relations ©^O, we use standard methods for slowly varying envelopes 
[21] to derive the dynamical equations for two pulses on a photon gas background, to 
obtain 



+ ck m ■ V )E 1 + 



> 01 



dt 1 ui ; 11 2k 

ck 02 ■V)E 2 



dt ^ J * 2k 



<02 



V 2 - (fcoi ■ V) 2 



V 2 - (fc 02 ■ V) 2 



Ex + Xckox [ 2 -<S g + ^a x S 2 ) E x = 0, (11) 
E, + Xrk (r2 ( 2 -£ g + X -a 2 gx j E, = 0.(12) 



I One notes that an overall factor two enters, because of the contributions from the terms where the 
wave field changes place with the complex conjugate term. However, this is compensated for by using 
the time averaged values of the fields rather than the peak amplitudes. 
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where the index refers to the non-interacting background state. We note that if the 
two pulses are co-propagating, their mutual interactions disappear, while for counter- 
propagating pulses we have a* = 4. 

Following Ref. [19] and using 0, we then derive the equations for the radiation 
gas on a background of two plane wave pulses (assuming that the rate of change in the 
direction of the pulses is slow) 



94 

dt 



c 2 v n 



4 



-A 




2 dt 



c 2 n • v Ki 




Yd-t +cn 



^ S 2 + c 2 (ft + £ 2 + f3^S^j 



V // 



, (13) 



and 



dn 1 



dt 



'-\s n 





vs. 



(14) 



where U is the momentum density of the photon gas. The coupled equations (jllj) - 
(H2J and (fT3j) - (fT4"j) describe the momentum exchange due to wave-wave scattering of 
the high frequency pulses and the low frequency gas photons. We note that Eqs. (HHJ)- 
(|14|) have been derived by taking moments of the Vlasov equation for the radiation gas 
photons. We have chosen to start from an effective Lagrangian for the photon-photon 
interaction, but one could in principle also start from the scattering cross section of 
the elastic photon-photon collisions, forming a collisional integral of the Boltzmann 
equation, and from there derive the fluid equations. However, for the purposes of this 
paper, the Vlasov approach is sufficient. 

Introducing S g = Sq+5S, where Sq (^> 5$) is the unperturbed radiation gas energy 
density in the absence of the pulses, and assuming that the unperturbed momentum 
density of the gas vanishes, we obtain from Eqs. (fT3 ]) -([T4" |) 



d 2 5<S 



df 



+ 



3 




^H- + c 2 V 2 
2\U X ) Ut 2 




+ c 2 |V^ 2 | 2 



1\J (%\<§2 



(15) 



Equation (fTo} together with Eqs. (fTT|) - (fr2|) form a Karpman-like system of equations 
[22-25]. 

In order to analyse modulational instabilities, we proceed along the lines of Ref. [26] 
and let 5$ = dcff exp[i(K ■ r — Qt)] and Sj = Sj expf^i^ • r — fit)} where j = 1,2. From 
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Eq. (|15[) we then obtain 
2 



5<f = -A4A ( S x + + /3 



(16) 



where A = (J? 2 + c 2 if 2 )/(-^ 2 + c 2 iT 2 /3). 

In the absence of the radiation gas, the interaction between co-propagating pulses 
vanishes. This is not the case when the radiation gas is present, as will now be 
shown. For parallel propagation, we have k 01 = k 02 = k , and V 2 — (k ■ V) 2 = Vj_. 
Furthermore, parallel pulses imply a\ = a 2 = 0, i.e., the Eqs. fTTf - fl"2j) take the form 



0_ 

Ft 



+ ck ■ V Ej + 



2 k, 



oj 



-V\E, 



l + -\e^{\E l \ 2 +\E 2 \ 2 + (3\E 1 \\E 2 \) 



Ea = 0, 



;i7) 



where we have used Eq. (fTB^) . The modulational instability of Eqs. (fTTj) will now be 



analysed, using the ansatz Ej 



'jjo + Ej)e iLp i l , where 



\EA < \E 



jo|> 



and where 



Ejo is a real constant vector. To lowest order, Eq. (fTTjl gives the frequency shift 
ifj = ~\ck j£a [l + |Ae A (E 2 Q + _Ef + j3E w E 20 )~\ , and to next order we obtain 

c 



9 t ^ 



2k, 



-V 2 : E< 



0.7 



+ El) ■ E w + (.E 2 + E 2 ) ■ E' 



20 



iSoQ — 0, 



;is) 



where the star denotes the complex conjugate, and Wj = (4/9)A 2 cfcoj^oA. 



To further simplify the analysis, we assume that Ej ■ E 



■jo 



EjEj 0} and that 



the polarisation of the perturbation vectors remain fixed to lowest order, and set 
&oi = ko2 = k . We let Ej = (Xj + iYj) exp[i(K ■ r — f2t)), and Fourier analyse 
Eq. (HH1) to obtain 



(Q-cK\\Y- 
where W 



cK\ 


'cK\ 


k 


_2k 


w x = w 2 , 4 



- ws n 



(n-cK\\) 2 + 



<fKl 



Fk n 



0. 



(19) 



e (E 2 + E 2 ), K \\ = k -K and K\ = K 2 - K 2 Thus, for 



Q = cK\\ + %T and F <C cK\\ , we find that 
1 



r 



-cK } 



16 K\ + 2K\ 
Y^K 2 -2K 2 



k 2 ' 



(20) 



where the dimensionless parameters i] = XS Q and i] p = XS p give the coupling strength 
for the radiation gas and the pulse, respectively. Thus, we see that the radiation gas 
acts as a mediator for the energy transfer between the two pulses, something that is not 
possible with the pulses alone. The positivity of T 2 requires that 

2 - 16k^T] r] p ) 



K , > Kl (3K 2 



> 0. 



■■ 2 (3K 2 ± + 16k 2 VoVp ) '- ~ (21) 

We furthermore notice that the growth rate diverges when K\\ approaches the "critical" 
value ±K±/y2, at which the approximation r <C cKu breaks down. On the other 
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hand, even if the parameters rj and r\ v are small, one can nevertheless obtain a large 
growth rate, provided that K\\ lies close to ±K±/y2. The phase speed at the critical 
wavenumber is c/v^3, i-e. the speed of sound in a thermal radiation gas. 

A stationary filamentation instability can likewise be obtained from Eq. (|17|). 
analogously to the modulational instability. If K\\ = ix, we obtain 

1 /l6 , , 

>t = 5^V"3**-#- (22) 

Notice that, contrary to the modulational instability, the filamentational instability 
growth rate does not have a singular behaviour at certain parameter values. 

The existence of a modulational and/or filamentational instability within a physical 
system indicates that the nonlinear evolution, either in time or in space, is nontrivial. 
When the time-dependence is significant, the modulational instability implies that 
any small perturbation will grow. We therefore expect the system to evolve towards 
a strongly inhomogeneous state. Naturally higher order nonlinearities are likely to 
saturate such a collapse scenario. On longer time-scales thermo-dynamical effects, 
outside the range of validity of our equations, will become important. This is likely to 
eventually destroy the inhomogeneous structures. Even if the system is stationary, the 
filamentational instability will single out a spatial direction along which inhomogeneities 
will grow. Here, a full nonlinear analysis requires numerical methods. We note that, 
in addition to the modulational and/or filamentational instabilities considered here, 
four-wave scattering involving only high frequency photons is possible. However, for 
dimensional reasons we expect the growth rates associated with such processes to be 
smaller than the ones found here, provided that rj > rj p . 

The situations in which photon-photon scattering can become important are indeed 
in the extreme range of what Nature may offer. Although the evolution of high power 
lasers has been rapid over the last decades, we are still not near the limits where the 
direct effect of lasers can give us important information of photon-photon scatterings. 
Rather, this has to be sought by more indirect means, such as the coherent interactions 
of cavity modes presented in Ref . [14] , or in astrophysical environments, such as close to 
magnetars (see, e.g., Ref. [27]), where field strengths close to the Schwinger limit may 
exist, and effects such as photon splitting or magnetic lensing, as found in Ref. [5-7,16], 
may take place. Moreover, the possibility of using recent high precision measurements of 
the cosmic microwave background, such as the Wilkinson Microwave Anisotropy Probe 
(see Ref. [28]), for the indirect detection of photon-photon scatterings has also been 
suggested in Ref. [19]. 

In this paper we have suggested a theory for pairs of electromagnetic pulses 
propagating on a radiation background. It has been shown that the dynamics of the 
photon pulses is governed by a set of coupled nonlinear Schrodinger equations, while 
the response of the background can be deduced from an acoustic wave equation driven 
by the pressure of the two photon pulses. The full system resembles a Karpman set 
of equations. Moreover, we have investigated the interaction between two parallel 
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propagating electromagnetic pulses via a background radiation gas. The gas acts 
as a mediator for the interaction, and the pulses exchange their energy through 
this background. In contrast to the case of no radiation background, when parallel 
propagation implies zero coupling, the radiation gas in our case gives rise to a system 
of equations possessing both modulational and filamentational instabilities. We have 
presented the temporal and spatial amplification rates for both instabilities, and have 
discussed the implications of these instabilities for the evolution of the system, both in 
laboratory settings and in astrophysical environments. 
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